Search for: All records

It is known that inhomogeneous second-order macroscopic traffic models can reproduce the phantom traffic jam phenomenon: whenever the sub-characteristic condition is violated, uniform traffic flow is unstable, and small perturbations grow into nonlinear traveling waves, called jamitons. In contrast, what is essentially unstudied is the question: which jamiton solutions are dynamically stable? To understand which stop-and-go traffic waves can arise through the dynamics of the model, this question is critical. This paper first presents a computational study demonstrating which types of jamitons do arise dynamically, and which do not. Then, a procedure is presented that characterizes the stability of jamitons. The study reveals that a critical component of this analysis is the proper treatment of the perturbations to the shocks, and of the neighborhood of the sonic points. 

Pressure Poisson equation (PPE) reformulations of the incompressible Navier–Stokes equations (NSE) replace the incompressibility constraint by a Poisson equation for the pressure and a suitable choice of boundary conditions. This yields a time-evolution equation for the velocity field only, with the pressure gradient acting as a nonlocal operator. Thus, numerical methods based on PPE reformulations are representatives of a class of methods that have no principal limitations in achieving high order. In this paper, it is studied to what extent high-order methods for the NSE can be obtained from a specific PPE reformulation with electric boundary conditions (EBC). To that end, implicit–explicit (IMEX) time-stepping is used to decouple the pressure solve from the velocity update, while avoiding a parabolic time-step restriction; and mixed finite elements are used in space, to capture the structure imposed by the EBC. Via numerical examples, it is demonstrated that the methodology can yield at least third order accuracy in space and time. 

Runge-Kutta time-stepping methods in general suffer from order reduction: the observed order of convergence may be less than the formal order when applied to certain stiff problems. Order reduction can be avoided by using methods with high stage order. However, diagonally-implicit Runge-Kutta (DIRK) schemes are limited to low stage order. In this paper we explore a weak stage order criterion, which for initial boundary value problems also serves to avoid order reduction, and which is compatible with a DIRK structure. We provide specific DIRK schemes of weak stage order up to 3, and demonstrate their performance in various examples.